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Abstract 

We have studied the precession mode, the rotational ex- 
citation built on the high-ZC' isomeric state, in comparison 
with the recently identified wobbling mode. The random- 
phase-approximation (RPA) formalism, which has been de- 
veloped for the nuclear wobbling motion, is invoked and the 
precession phonon is obtained by the non-collective axially- 
symmetric limit of the formalism. The excitation energies 
and the electromagnetic properties of the precession bands in 
^^**W are calculated, and it is found that the results of RPA 
calculations well correspond to those of the rotor model; the 
correspondence can be understood by an adiabatic approxi- 
mation to the RPA phonon. As a by-product, it is also found 
that the problem of too small out-of-band B{E2) in our pre- 
vious RPA wobbling calculations can be solved by a suitable 
choice of the triaxial deformation which corresponds to the 
one used in the rotor model. 

1. Introduction 

In this talk, we present a recent progress of our study 
on the wobbling and precession modes in nuclei. These 
keywords, wobbling and precession, represents the mo- 
tions of classical tops. They are quite interesting be- 
cause they describe three-dimensional (3-D) rotational 
motion, and so related to a fundamental question: How 
does a nucleus rotate as a 3-D quantum object? 

We do not know how strictly these two are distin- 
guished, but in the following we use wobbling for mo- 
tions of triaxial-body and precession for those of axially- 
symmetric body. In the lowest energy motion, the top 
rotates about one of the principal axes with largest mo- 
ment of inertia, but when excited the angular momen- 
tum vector tilts from this axis in the body-fixed (intrin- 
sic) frame. Then, looking it from the laboratory frame, 
the body wobbles or precesses, and that is why these 
names came from. In the classical mechanics, these two 
are similar; actually the difference is that the trajectory 
of the angular momentum vector is a pure circle in the 
case of precession, while it is an ellipse in the case of 
wobbling. 

However, the atomic nucleus is a quantum system and 
the situation is dramatically changed: The collective 
rotation cannot occur about the symmetry axis. Thus, 
the quantum spectra corresponding to the wobbling and 
precession motions are completely different. In the case 
of wobbling, the spectra associated with an intrinsic con- 



figuration are composed of multiple rotational bands; 
the lowest (yrast) one represents the uniform rotation 
about the main rotation axis, the first (onc-phonon) 
excited band represents a quantized motion of tilting 
angular momentum vector, and so on (more excited 
band with more tilting angle). These {AI = 2) rota- 
tional bands corresponds to a collective rotation about 
the main rotation axis with largest moment of inertia. 
Since the excitation of phonons, or tilting the angular 
momentum vector is another type of rotation about the 
axis perpendicular to the main rotation axis, these ex- 
citations form [AI = 1) vertical sequences. In this way, 
the wobbling motions show a complicated band struc- 
ture, the horizontal rotational sequences and vertical 
phonon-like excitations. 

On the other hand, the angular momentum along the 
symmetry axis (the main rotation axis) is generated by 
quasi-particle alignments in the high-if isomeric con- 
figurations, and no collective rotation exists about it 
in the case of precession. There are no horizontal se- 
quences leaving only one vertical band for each intrinsic 
configuration, which is nothing but a collective rotation 
about the perpendicular axis to the high-i^ angular mo- 
mentum. 

Recently, the wobbling phonon spectra have been 
identified among the triaxial superdeformed (TSD) 
bands in Lu nuclei T, 'S], up to two-phonon excita- 
tions. By using the microscopic framework, the random- 
phase approximation (RPA), we have studied the wob- 
bling phonon in the Hf-Lu region |3, Q . The precession 
bands are rotational bands excited on the prolate high- 
K isomers and have been known for many years, see 
e.g. We have recently investigated the properties of 
the precession modes applying the same RPA formalism 
in comparison with the wobbling modes (5|, |^ . The con- 
tent of this talk is largely based on the results of , and 
a further development on the out-of-band B{E2) transi- 
tion probability of wobbling phonon excitation (see §3). 

2. Precession mode as a phonon 

The study of wobbling and precession modes has a 
long history (see references quoted in 6]). Our recent 
works on the wobbling motions rely on the microscopic 
framework developed in 1979 by Marshalek jgj. In al- 
most the same time, a very important work for the pre- 
cession mode have been done here in Lund in 1981 by 
Andersson et al. (c.f. also Refs. [13 [ill). In these 
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works, was used the microscopic RPA formalism, which 
is known to be suitable to describe vibrational excita- 
tions. The reason why the RPA is employed for such 
specific rotational motions as the wobbling and preces- 
sion modes can be easily understood in the case of pre- 
cession band, i.e. the high- if rotational band. 
The spectra of high-ii' rotor is well-known: 
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Putting n = I — K and assuming K is large, it is easy to 
see that this leads to an approximately harmonic spec- 
trum, 
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where Wprec = K/J'± is the precession phonon energy. 
The remaining anharmonic term is of order (l/K) and 
can be neglected at the high-X limit. This harmonic 
picture is also valid for the E2 transitions, 



B{E2) = (5/167r) Ql {IfK2{)\UKf 
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By taking the high-iiT asymptotic limit in the Clebsch- 
Gordan coefficient. 
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in which the one- phonon transition is order (l/K) and 
the direct transition to two-phonon state is {1/K'^), so 
that the direct two-phonon transitions are prohibited in 
the high-JiT limit. 

Now it is clear that the phonon picture is valid and 
the RPA formalism can be used to describe the pre- 
cession mode microscopically, which has been done by 
Andersson et al. ||9|. The key of their work is to adopt 
the so-called symmetry-restoring separable type interac- 
tion: The residual interaction used in the RPA is con- 
structed in such a way that the rotational symmetry 
broken by the deformed mean-field hamiltonian is recov- 
ered. This uniquely determines the complete form of the 
residual interaction, and there is no adjustable param- 
eter in the RPA calculation. The dispersion equation 
is simply given as lo S{uj) = 0, where the uj ~ solu- 
tion is the rotational symmetry-restoring mode, i.e. the 
Nambu-Goldstone (NG) mode, J± — Jy ± iJz- Here 
and in the following we assume the x-axis as the high- 
ly alignment axis. The equation for non NG modes can 
be formally cast into the same form as Wproc = K/J^±, 
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but with microscopically defined moment of inertia. 
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However, the inertia J^^\ujprec) is energy-dependent 
and actually Eq. |0J is a non-linear equation to obtain 
the RPA eigen-energy. Electromagnetic transitions can 
be calculated as the squared amplitudes of the RPA 
phonon X^,.^^ with respect to appropriate transition op- 
erators Qx^: 
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Here it should be stressed that this RPA formalism of 
Andersson et al. is obtained by taking the non-collective 
axially symmetric limit (7 = 60° or —120° in the Lund 
convention) of the RPA wobbling formalism of Mar- 
shalek 6] . In this way, the precession and wobbling can 
be considered as similar kinds of collective excitation 
modes. 
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FIG. 1: Excitation energies of the one-phonon precession 
modes excited on high-iC' configurations. Calculated ones 
are denoted by filled circles connected by solid lines, and 
experimental ones by crosses. Taken from Ref. |6|. 



The results of RPA calculations for precession modes 
in I'^^W 7| are presented in Fig. ^ In this nucleus 
there observed many high-iC isomers, on which the ro- 
tational bands exist: Eleven isomers investigated are 
ranging from four-quasiparticle states (2i/-27r) to ten- 
quasiparticle states {Qv-Att). We have used the Nilsson 
potential as a mean field with appropriate deformation 
parameters and pairing gap parameters. Except for four 
cases, = 18~, 25+, 28~, 29+, the one-phonon exci- 
tation energies are well reproduced. As for these four 
isomers, the precession energies are too small; namely, 
the calculated moments of inertia are too large. We 
found that the overestimation of the calculated moments 
of inertia is due to the proton contributions in these 
four configurations, which include the 7r[541]l/2~ Nils- 
son state originated from the 77/19/2 high-j decoupled 
orbit. Apparently the contribution from this orbit is 
overestimated in the present calculations, see Ref. 
for more detailed discussions. 

As for the electromagnetic transitions, one can di- 
rectly compare the calculated and measured transition 
probabilities, but here we compare them in a different 
way. As is well-known, the rotor model gives simple for- 
mula for transition probabilities, e.g. Eq. © for B{E2). 
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FIG. 2: {qk — (7H)-factors for high-Tf precession bands. 
Those calculated by the RPA are denoted by filled circles 
connected by solid lines, while those by the mean-field ap- 
proximation by filled triangles connected by dotted lines. 
Those extracted from the experimental data are shown by 
crosses. Taken from Ref. 0. 



Thus, by using the asymptotic forms of the Clebsch- 
Gordan coefficients, B{E2) and B{Ml) are parameter- 
ized in the high-/^ limit as. 
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On the other hand, they are calculated by Eq. Q in the 
RPA formalism. By equating these two expressions, we 
define Q moment or gK — gn factor calculated within 
the RPA. Namely, we parameterize the results of RPA 
calculations in the same way as the rotor model, and 
they are compared with experimentally extracted ones 
in Fig. 121 for B{M\). In this figure we included the 
results of the simple mean-field approximation, namely 
using gK calculated as {fix)/{Jx) (a*^ is the Ml opera- 
tor) with each high- if state, and a common value of ga 
calculated in the same way but with the ground state. 
As it is clearly seen in Fig. [3 the RPA calculation gives 
a much better description of B{AI1) values. As for the 
E2 transitions, there is no experimental data available, 
but by comparing the calculated RPA Q moments with 
the mean-field estimate of Q moments it is found that 
they well coincide except four configurations, which in- 
clude the 7r/ig/2 decoupled orbit (the figure is not shown, 
see Ref. Q). 

The results for the excitation energies and transition 
probabilities indicate that the high- if rotor model pic- 
ture is realized if the model parameters, moment of in- 
ertia, Q moment and g-factors, are calculated appro- 
priately by means of the RPA. This fact can be natu- 
rally understood by using an adiabatic approximation 
studied in Refs. ^ [ill; i.e. in this approximation, 
(J+ + KiO+)B,PA/V2K, while the NG mode 
= J_ / \J2K. By using this form it is easy to 
confirm, for example, ((5o)mean-fieid ~ ((9o)rpa; the sit- 
uation is more complicated for the g-factors, see Ref. 
for more detailed discussions. 
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3. Out-of-band transition of wobbling phonon 

Now let us come back to the case of wobbling mode. 
It has been shown that the RPA calculation well corre- 
sponds to the rotor model picture in the case of preces- 
sion, especially for B{E2). It can be shown that a simi- 
lar consideration in term of the adiabatic approximation 
can be also applied for the wobbling mode. However, the 
B{E2) value of our RPA wobbling calculations in Lu nu- 
clei were too small, by about 1/2 — 1/3, compared with 
the experimental data, whose magnitudes are well re- 
produced by the rotor model. This has been a serious 
problem for us, see Fig. Eland Refs. 0,11. Considering 
the results of the precession modes, however, it is dif- 
ficult to understand the discrepancy between results of 
the RPA calculation and the rotor model. Therefore, 
we looked for the reason why the RPA calculation gave 
such smaller values. 
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FIG. 3: Upper panel: The relation between 7(Nils) and 
7(dens) for the TSD band in ^^'^Lu. Lower panel: Compar- 
ison of the out-of-band to in-band B{E2) ratios; symbols 
are the experimental data, three curves are corresponding 
results of RPA calculations, in which the dashed line is the 
result of Ref. H, see text for detail. 



We have found that there are two reasons: One is re- 
lated to the model space of the RPA calculations, and 
the other is to the definition of triaxiality parameter 7 in 
the Nilsson potential. As for the model space, we used 
the 5-major shells but it was not enough; this effect, 
however, is only about 20% and not the major effect. 
More important is the value of the 7 parameter. It is 
believed that the 7 value is about 20° in the TSD bands 
in the Lu region. We have used 20° for the 7 in the Nils- 
son potential, 7(Nils) = 20°. It is found, however, that 
the 7- value defined by the density distribution, 7(dens), 
is only about 10°. Here, 7(dens) is defined by the ex- 
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pectation values of the quadrupole operators (the Lund 
convention for sign of 7), — tan~'^((522^^)/(Q20^'')' i-^- 



7 (dens) = tan 



-1 a=l 



(8) 



($:(2.^-x^-2/2).) 
0=1 



Of course, the triaxiahty of the rotor model should be 
that of the density distribution. By using an assump- 
tion for moments of inertia, J^y = J7z, in the triaxial 
rotor model, the out-of-band to in-band B{E2) ratio 
depends on the 7 parameter like B{E2)out/ B{E2)in oc 
tan^ (7 + 30°) Even in general cases with Jy ^ Jzi 
this 7-dependence approximately follows. Then it is 
easy to see that the difference of this ratio at 7 = 10° 
and 20° is about factor two. The 7 value corresponding 
to 7 (dens) = 20° in the density distribution is about 
7(Nils) « 30° in the Nilsson potential (see Fig.EJ). The 
results of the RPA calculation using different 7-values 
are shown in Fig.|3| where the relation between 7 (dens) 
and 7 (Nils) is also depicted in the upper panel. The 
dashed line with the model space of the 5-major shells 
is the result of Ref. As is shown in the figure, if one 
uses the full model space and 7(Nils) = 30°, the B{E2) 
ratio comes up to the correct magnitude just like in the 
case of the rotor model |13]. 

The discrepancy between 7(Nils) and 7(dens) is just a 
matter of definition, see e.g. Appendix of Ref. M^. For 
the pure harmonic oscillator potential, the definition of 
7 (Nils) gives 



7 (Nils) = tan 



-1 \/i{uoy ~ i^x) 
2lo, — [jJx — uj,, 



(9) 



The selfconsistency of the potential (Mottelson condi- 
tion), (I]a=i(2;fc)a) = 2;,?/,z), is approxi- 
mately satisfied in the Nilsson potential, and then by 
Eq. 0, 



7(dens) « 7(self) = tan" 



V3(l/c.g-l/c.g) 

2/0.2 - _ 1/^2 
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where 7(self) is the triaxiality based on the shape (el- 
lipsoid) of an equi-potential surface for the anisotropic 
harmonic oscillator. In this way the relation between 
7(Nils) and 7(dcns) in the upper panel of Fig. Ocan be 
naturally understood. 

4. Summary 

The RPA precession formalism by Andersson et al. 
can be obtained by a non-collective axially-symmctric 
limit of the RPA wobbling formalism by Marshalek. By 
using this formalism, it is shown that the RPA calcu- 
lation gives fairly good agreements for both the excita- 
tion energies and B{MX) transitions for the precession 
phonons on high-iC isomers in ^^^W. This result can be 
interpreted by an adiabatic approximation, and shows a 
good correspondence between the RPA calculation and 
the rotor model; especially B{E2) or Q moments, and 
B{M1) or 5-factors 0. 

There is an important feed back to the calculation 
of recently observed nuclear wobbling motions. The 
problem of small B{E2) ratios in our previous RPA 
calculations can be solved if one uses a proper value 
of the triaxiality parameter 7(dens) w 20° in the den- 
sity distribution. This does not completely solve the 
problem, however, because the Nilsson-Strutinsky cal- 
culation gives minima at 7(Nils) « 20° in the Nilsson 
potential, which corresponds to 7(dens) ~ 10°. Thus, 
it raises a more fundamental question; why the Nilsson- 
Strutinsky calculation does not provide enough triaxial- 
ity, which is required for explaining the measured B{E2) 
ratio of the wobbling phonon band. Finally, it is shown 
that the RPA can describe the 3-D rotational motion 
(in the small amplitude approximation) at least as the 
same level as the rotor model. In other words, the rotor 
model may be justified by the microscopic RPA calcu- 
lations. 
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